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A constraint on the viable f{R) model is investigated by confronting theoretical predictions with 
the multipole power spectrum of the luminous red galaxy sample of the Sloan Digital Sky survey 
data release 7. We obtain a constraint on the Compton wavelength parameter of the f{R) model on 
the scales of cosmological large-scale structure. A prospect of constraining the Compton wavelength 
parameter with a future redshift survey is also investigated. The usefulness of the redshift-space 
' distortion for testing the gravity theory on cosmological scales is demonstrated. 
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I. INTRODUCTION 

Experimental tests of gravity on the scale of the Solar System show good agreement with predictions of general 
relativity (e.g., [ij). The nature of the Newtonian gravity is the attractive force, which naturally predicts a decelerated 
expansion of the universe. Contrary to this expectation, it has been discovered that our universe is undergoing an 
f~| ' accelerated expansion epoch 0-13l- Though the accelerated expansion is explained by introducing a cosmological 
CIh" constant, its small but nonzero value cannot be explained naturally 0]. The problem might be deeply rooted in the 
Q ' nature of fundamental physics. 

5^ . This problem has attracted many researchers, and many works have been done, both theoretically and observation- 
' ally. As a generalisation of the cosmological constant, a dynamical field, called the dark energy model and its variants, 
Ci . are proposed to explain the accelerated expansion of the universe (see |6| and references therein). As an alternative 
to the dark energy model, modification of gravity may explain the accelerated expansion. General relativity is not 
considered to be the complete theory, because its quantum theory cannot be formulated in a well defined manner. 
The theory of gravity might need to be reformulated within a more general framework. 

From the observational point of view, the constraint on the gravity theory on cosmological scales has not been well 
CO . investigated, compared with the constraint on the scales of the Solar System. Many future projects to produce large 
■ galaxy surveys are in progress or planned which aim to explore the nature of the dark energy. These surveys 

I are useful for testing the theory of gravity at cosmological scales (e.g., [l^l)- The dynamical dark energy models may 
have similar expansion rates as models of modified gravity, but predict different histories for the growth of structures. 
, The key to testing the gravity theory is the measurement of the evolution of cosmological perturbations, as many 
authors have concluded recently |14| - [25j . 

The cosmic microwave background anisotropics are useful for investigating the cosmological perturbations through 



the measurements of the integrated Sachs- Wolfe effect or the lensing effect on the angular power spectrum 
Imaging surveys of galaxies are also useful through the weak lensing statistics or cluster number counts (27l . 



261. 



Similarly, redshift surveys of galaxies are helpful for testing gravity [29|-|3J|. In the present paper, we revisit the 
' problem of testing the gravity theory through a measurement of the multipole power spectra in the sloan digital sky 



survey (SDSS) luminous red galaxy (LRG) sample [3l|. Measuring the multipole power spectra is a wa y to estimate 
the redshift-space distortions, which reflects the linear growth rate of the matter density perturbations [35l- l37| . 

Many authors have investigated the clustering nature of the SDSS LRG sample [38l - l47| . In the references [48l [iot , 
recent results on LRGs from the SDSS data release (DR) 7 are reported. In the reference [i^, a test of gravity is 
considered using the observed anisotropic correlation function. Three of the authors of the present paper have shown 
that the SDSS LRG sample is useful to test the gravity theory by measuring the quadrupole power spectrum of galaxy 
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distribution, which represents the redshift-space distortions 'si'l. In the present paper, we revisit the issue of testing 
the gravity theories on the cosmological scales using the SDSS LRG sample of the DR 7, especially focusing on the 
f{R) gravity model. 

The f{R) models proposed in 'so'-'SS^ are viable models of modified gravity, which include some function of the 
Ricci scalar, f{R), added to the Einstein Hilbert action. As the modification of gravity involves the introduction 
of extra degree of freedom in general, one must be careful with the resulting behaviour. Furthermore, any theory 
must reduce to the general relativity on the scales of the Solar System. In the f{R) model, the general relativity is 
supposed to be recovered by the chameleon mechanism [58, 59], which hides the field of the extra degree of freedom 
because the mass of the field becomes large for a dense region. The cosmological bounds on the f{R) model have 
been investigated with the cosmic microwave background anisotropics (60j and also using the abundance of galaxy 
clusters [6l|. However, our approach is based on the redshift-space distortion [ooj . 

This paper is organised as follows: In section 2, we briefiy review the f{R) model and the characteristic evolution 
of the matter density perturbation. In section 3, we present our results for the multipole power spectrum of the SDSS 
LRG sample of the DR 7. In section 4, cosmological constraint is discussed by confronting the observed multipole 
spectra with the theoretical predictions. In section 5, a prospect of constraining the f{R) model is discussed on 
the basis of the Fisher matrix analysis, assuming a future large redshift survey. Section 6 is devoted to summary 
and conclusions. Throughout this paper, we use units in which the velocity of light equals 1, and adopt the Hubble 
parameter Hq = 100/ikm/s/Mpc with h — 0.7. 



II. f{R) GRAVITY MODEL 

In this section, we briefly review the f{R) model, proposed in the references [50l - [53| . In general, higher derivative 
terms are expected in the low energy effective action of gravity. Inspired by this, the f{R) model introduces some 
function of the Ricci scalar f{R), adding to the Einstein Hilbert action. We consider the theory defined by 

S = I d^x^g {R + f(R)) + Sm, (1) 

where S'm is the action of the matter. Many aspects of the f{R) model have been investigated; see e.g. [13, [IHl 
for a review (cf., [sgI ISTj). We assume that the chameleon mechanism is responsible for the recovery of the general 
relativity on the Solar- System scales. The chameleon mechanism is a nonlinear effect. Recently, the effect on the 
quasi-nonlinear power spectrum is investigated based on the perturbative approach or the numerical simulations (62| - 
1651 . This nonlinear chameleon effect becomes influential in the nonlinear regime. In the present paper, however, we 
can neglect the nonlinear chameleon effect because we need to consider only rather large scales, k ^ 0.2/iMpc~^. 

For the viable model, the function f{R) must satisfy some conditions. We consider the model where the asymptotic 
form of f{R) can be expressed by 



f{R) -2A 



R, ^ 



(2) 



where A is the cosmological constant, n is a constant that specifies the f{R) model, and Rc is also a constant with 
the same dimension as that of the Ricci scalar. The background expansion of this f{R) model is well approximated 
by that of the ACDM model. 

It is known that the additional term f{R) involves the introduction of an extra degree of freedom. Namely, 
fji = df /dR corresponds to the extra degree of freedom, which behaves like a scalar field. From the above action, one 
can derive the equation for fj^, 

l{R + 2f-Rf,)+^-f- 

where p and P are the energy density and the pressure of the matter, respectively. If we regard the right hand side 
of equation ^ as the derivative of the effective potential, dVcs/df]^ , the mass of fn can be read 



V^V/fl = -{R + 2f- Rfn) + —hp + 3P), (3) 



2 d^^eff Ifl + fR 



The viable f{R) theory satisfies f R, and |/i^| <C 1. Assuming R/rr ^ Ij the mass of the extra degree of freedom 
is 

- (5) 
3 JRR 
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FIG. 1: (a, Left) Di{a)/a as a function of the scale factor. The soUd curve is the ACDM model with f2o = 0.28. The dashed 
curves are for the f{R) model with the wavenumbers fc/(/iMpc^^) — 0.2, 0.1, 0.05 from the top to bottom, respectively. Here 
we adopted the model n = 1 and kc = 0.05/iMpc~^. (b. Right) Same as (a) but for the growth factor / — d\nDi/d\na. 



where ffm ~ <P f /dE? . Thus, fB,R > is required to avoid the extra degree of freedom to become tachyonic. This 
extra degree of freedom mediates an attractive force, and modifies the gravity from the range determined by the 
Compton wavelength A 1/m. From Eq. ([2]), we have 

In the subhorizon limit, the matter density perturbation follows (e.g., (66j and references therein), 

J + 2-(j-47rGoff(a,fc)p(5 = 0, (7) 
a 



where 



G'off(Q,fc) _ 1 fcVa^ 

G ^ 3fc2/a2 + l/(3/iii?)' 



and the dot denotes the differentiation with respect to the cosmic time. 
Instead of Rc, we introduce the parameter by 

17o/a^ + 4(1 - f7o) 



n f^o + 4(1 - fio) 



(9) 



where represents the wavenumber corresponding to the Compton wavelength at the present epoch. Thus, the f{R) 
model is specified by n and k^.. The growth factor can be obtained by solving Eq. ([7]), which we denote by Di{a, k). 
The growth rate is given by / = (ilnZ?i(a, k)/ d\na. 

In the Einstein de Sitter background universe, the evolution of the density perturbation can be solved analytically 
[68l | . Two of the authors of the present paper investigated characteristic features of the evolution of the growth rate of 
the /(i?) model, both numerically and analytically in the reference [67| . In the present paper, we solve the evolution 
equation ([7]) numerically (cf. [gl, [t^]). Figure [T] shows the growth factor divided by the scale factor (left) and the 
growth rate (right), respectively, as a function of the scale factor. The solid curve is the ACDM model with the density 
parameter rip = 0.28. The dashed curves are for the f{R) model with different wavenumbers fc/(/iMpc~^) = 0.2, 0.1, 
and 0.05, respectively. Here the /(i?) model assumes n — 1 and kc — O.OSft.Mpc""'^. Due to the modification of the 
gravity the growth factor and the growth rate are enhanced, and this enhancement is scale-dependent. 
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FIG. 2: The mean number density of galaxies, n, as a function of the redshift z of the SDSS LRG sample, where we adopted 
the ACDM model with ilo = 0.28 for the distance-redshift relation s = s\z\. 



III. MULTIPOLE SPECTRUM OF THE SDSS LRG SAMPLE 



The multipole power spectrum Pi{k) is defined by the coefficient of the multipole expansion of the anisotropic 
power spectrum P{k, /i), 



Pik,f,)^ J2 Piik)Cei^x)i2£+l), 



(10) 



where CeXn) are the Legendre polynomials, /i(= cos 9) is the directional cosine between the line of sight direction and 
the wavenumber vector k. Note that our definition of the multipole spectrum Pe{k) is different from the conventional 
one by the factor 21 + 1 [H, IH, [7l| . Here the Legendre polynomials satisfy the normalisation condition. 



+1 



21+1 



(11) 



The monopole Po{k) represents the angular averaged power spectrum, which is what we usually mean by the power 
spectrum; the quadrupole P2(k) represents the leading anisotropy in the power spectrum due to the redshift-space 
distortion. The hexadecapole P^ik) represents a different aspect of the redshift-space distortion. In the present paper, 
we focus on the monopole and quadrupole spectra. The quadrupole spectrum reflects the peculiar velocities of the 
galaxies [3^[3^[7i|]. Those peculiar motions can be used to test the gravity theory on cosmological scales. 

Pioneering works on the measurement of the quadrupole spectrum was carried out by Cole, Fisher, and Weinberg 
[35} and Hamilton [36| using the IRAS galaxy survey catalogue. Cole et al. presented a systematic method to estimate 
the quadrupole power spectrum through the anisotropic power spectrum [35| . The method was applied to the Two 
Degree Field (2dF) galaxy survey to estimate the /3 factor. Hamilton obtained the quadrupole power spectrum by a 
transformation of the correlation functions [s^. In the present work, however, we adopt a different method to estimate 
the quadrupole power spectrum [75| . Our method is in line to the widely used way to estimate the monopole power 
spectrum j72l [73| , and allows us to obtain the multipoles of the redshift-space power spectrum without evaluating 
the correlation function or the anisotropic power spectrum. In Ref. ^l|, we applied the method to the SDSS LRG 
sample from DR 6 to test the general relativity on cosmological scales. In the present paper, we revisit this problem 
with the SDSS LRG sample of DR 7 ^T^. 

Our LRG sample is restricted to the redshift range z — 0.16 — 0.47. In order to reduce the sidelobes of the 
survey window we remove some noncontiguous parts of the sample (e.g. three southern slices), which leads us to 
~ 7150 deg^(= A^) sky coverage with a total of TV = 100157 LRGs. The data reduction procedure is the same 
as that described in [3^. In this power spectrum analysis, we adopted the spatially flat Lambda cold dark matter 
(ACDM) model distance-redshift relation s = s[z], which is consistently chosen when comparing with theoretical 
prediction. 

The strategy to measure the multipole power spectrum is the same as that described in [tH] • We adopt the estimator 
of the multipole power spectrum for the discrete density fleld of the galaxy catalogue, as follows, 



Pt{k) 



AVfe 



d3fc(i?,(k)-^,(k)). 



(12) 
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FIG. 3: Po{k) of the SDSS LRG sample, where we adopted the distance-redshift relation s — s[z] of the ACDM model with 
fio ~ 0.28. The dark (black) points correspond to the DR7, while the light (green) ones to the DR6. The dashed and dotted 
curves show the f{R) model with n — 1/2, adopting a scale-dependent bias (case 1 in Eq. (|23|) ) and ctv — 350km/s. The dashed 
curve is for kc = l/iMpc~^, while the dotted curve is kc = 10~^/iMpc~^. The cosmological parameters are fio = 0.28, h — 0.7, 
and Us — 0.96 (primordial spectral index), and the amplitude of the perturbation is determined so as to be ug = 0.8 in the 
limit of infinitely large kc- 
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FIG. 4: P2{k)/ Po{k) of the SDSS LRG sample. The meaning of the points corresponds to those of Fig.[3l The dashed (dotted) 
curve is the theoretical prediction of the f{R) model with n — 1/2, kc — l/iMpc~^ (10~''/iMpc~^). The parameters of the 
bias model and CTv are the same as those of Fig. [3] The other cosmological parameters and the amplitude of the primordial 
perturbation of the f{R) model are also the same as those of Fig. O 
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where AVfc is the shell in the Fourier space and 



N 



A^rnd 



h ji 
N 

5,(k) =A-i(l + a)^7/;(s,,,k)£,(s,, • k), 



(13) 
(14) 



where s^j^ (sji) is the position of galaxies (random sample), ^ is the weight factor, which we take ■(/' = l,/i = s- kis 
the directional cosine between s(= s/|s|) and k(= k/|k|), a = N/N^-nd in our case is 0.05, and A is determined by 



A = 



(isn^(z)V'^(s, k). 



(15) 



Here the integral in the expression for A means the integration over the whole survey volume, and n{z) is the mean 
(comoving) number density of the galaxies. The error for the estimator Pi{k) is given by the variance [75|, 



(AP,(fc)2) 



with 



QKk) 



AVk 



dkA- 



AV,. 



s{Zn 



dsn'*(z)V/(s, k)[P(k, s) + l/n{s)YCi{s ■ k) 



(16) 



(17) 



Here we have assumed a <C 1. The covariance between the errors of different multipole spectra {APe{k)AP£i{k)) can 
be evaluated with the same formulae (fT6|) and (|17p . but only replacing £^(8 • k) by £f(s • k)Ce'{s ■ k) in (|17p [9lj |. In 
our analysis we adopt ip{s,'k) = 1. Figure [2] shows the mean number density as a function of z, when assuming the 
ACDM with Hq — 0.28 for the distance-redshift relation s = s[z]. 

Figure [3] compares the observed monopole power spectrum and our theoretical model. The dark(black) points with 
error bars in figure El show the monopole power spectrum of the DR7. The light(green) points are the previous results 
for the DR6 [31|. The dashed and the dotted curves represent the f{R) model with n = 1/2, with the scale-dependent 
bias model of case 1 (see the next section for details). The dashed curve is for kc = l/iMpc^^, while the dotted one 
for 10~'^/iMpc~"'^. The cosmological parameters are ilg = 0.28, h — 0.7, rig = 0.96 (primordial spectral index). The 
amplitude of the primordial perturbation is chosen to be erg = 0.8 in the limit of infinitely large kc- The Smith's 
nonlinear fitting formula [s^ is adopted. One can see that Po{k) can be fitted with our theoretical model, by choosing 
suitable bias parameters. 

Figure m plots P2{k)/ PQ{k). The meaning of the points and the parameters of the curves corresponds to those of 
Fig. [31 This figure shows that the quadrupole power spectrum can be used to constrain the f{R) model. Also it is 
clear that the long Compton wavelength model doesn't fit the data. 



IV. COSMOLOGICAL CONSTRAINT 

In order to investigate the cosmological constraint on the f{R) model from the multipole spectra, our theoretical 
model needs to include nonlinear effects. In the present paper, for simplicity, we adopt the following model for the 
galaxy power spectrum [t^ [t^] , 

Pgai {k, fi, z) = {b + fn^)^Pni{k, z)V KM , (18) 

where Pn\{k, z) denotes a nonlinear matter power spectrum, Vlkfiav] is the damping factor due to the Finger of God 
effect, and cr^ is the pairwise velocity dispersion. Assuming an exponential distribution function for the pairwise 
velocity, e~^l"" 1/*^^ / y^gy, where Vi2 is the pairwise peculiar velocity projected along the separation of a pair, the 
damping function is [13 (cf. Hlzil), 



T> [tTvfc/Lt] 



l+?2fc2^2/2 



(19) 
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FIG. 5: Ax^ on the kc — (Jv plane. Here we adopted the model n = 1/2. The other parameters are fio ~ 0.28, h — 0.7, 
Us = 0.96. The normalisation is fixed (Tg = 0.8 in the limit of large kc- The Peacock and Dodds's nonlinear fitting formula is 
used for the thin curves, while the Smith formula is used for the thick curves. Solid (dashed) contours correspond to /S.x^ — 6.2 
(2.3). The left panel adopted Eq. H25|) . while the right panel the covariance matrix from the mock catalogues. 
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FIG. 6: The correlation matrix, Eq. ((26}, for £ = (left) and 1 = 2 (right), respectively, from 1000 mock catalogues. 



with (Jv (^v/Hq. In this case, we have 
1 
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FIG. 7: Ax'^ on the kc — n plane, which we evaluated with Eqs. (1241) and (|25l) . For each pair of kc and n, the minimum value 
of is computed by fitting the bias parameter and Uv Other parameters are fixed f^o ~ 0.28, h — 0.7, and Us — 0.96. The 
normalisation of the primordial perturbation is chosen so as to be as = 0.8 (a) as = 0.82 (b), and as = 0.78 (c), in the limit 
of large fee. The panels (a)-(c) adopt the bias model of case 1. The panel (d) is the same as (a) but with bias model of case 
2. Solid (dotted) contours correspond to Ax^ = 6.2 (2.3). Almost overlapping thin and thick curves assume the Peacock and 
Dodds's formula and the Smith's formula, respectively. 



from Eqs. (|18l) and (|19p . For the nonlinear matter power spectrum, Pn\{k, z), we adopt the fitting formulas by Peacock 
and Dodds (Slj or by Smith et al. [8^. For the bias, we consider the following scale-dependent forms, 

(case 1) 




b{k)={ V-.,..v^p. ( / . x2 , (23) 



0.1/iMpc"\ 

1 I r I + ^2 ( T ) (case 2) 

VO.l/iMpc^V Vo.l/iMpc^V 

where 6o, 6i, 62, and a are the fitting parameters. 

Our strategy is the following. We use the monopole and quadrupole spectra in the wavenumber range 0.02/iMpc"^ < 
ki < 0.2/iMpc~^, and compute the chi squared 

= E E(^^(^*) - P^'''ik^))Cu}{k^, k,){P,, (fc,) - P^}^^{k,)), (24) 

where Pf"^{ki) is the observed power spectrum and Cu'{ki,kj) ~ {APi{ki)APi'{kj)) is the covariance matrix. Here 
the covariance of the errors of the monopole and quadrupole spectra is taken into account, however it does not affect 
our results quantitatively. 

The left panel of Fig. [5] shows the contours of Ax^ on the kc — CTv plane, where we used the covariance matrix from 
section 3, 

Cu,{h,k,) = {APeih)APf,ik,)) S^,. (25) 

The one-sigma (dashed curve) and two-sigma (solid curve) contour-levels are given, respectively. Here the chi squared 
is computed to minimise (1241) by fitting the bias parameters bo, 61, 62, or a, for each value of kc and cr^. The other 
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FIG. 8: The same as the Fig. [7] but with the covariance matrix from the mock catalogues. 
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FIG. 9: The same as the Fig.[7]but with the covariance matrix from the mock catalogues and the redshift-space power spectrum 
(|27|) . Only the curves with Peacock Dodds's formula for the nonlinear matter power spectrum are plotted. 
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parameters are fixed n = 1/2, fJo = 0.28, J7b = 0.044, Ug = 0.96, and h = 0.7. For P„i(fc, z), we adopted the Peacock 
and Dodds's formula [8l| (thin curve) and the Smith formula 182)] (thick curve), respectively. The redshift is fixed 
to z = 0.3, which is typical for the LRG sample. The amplitude of the matter power spectrum is fixed so as to be 
as = 0.8 in the limit of infinitely large fee, i-e., in the limit of the ACDM model. 

For comparison, the right panel of Fig. [5] shows the contours of Ax^, which take the correlation of the errors of 
different wavenumbers into account by evaluating Eq. ([24]), with the covariance matrix obtained from mock catalogues. 
Due to the inclusion of the correlation of errors of different wavenumbers, the constraint becomes weaker compared 
with the left panel. 

In the right panel of Fig. [SJ we obtain the covariance matrix by using mock catalogues, which were built by following 
the procedure described in the reference [s^. First, we generate density field using a second order Lagrangian 
perturbation calculation. Then, we perform Poisson sampling of the generated density field so as to end up with 
a galaxy sample that has a clustering strength enhanced by a bias and a number density equal to the observed 
LRG sample density. We then extract the catalogue by applying the radial and angular selection function. We have 
checked that the mock catalogues have the amplitude of the monopole and quadrupole power spectra consistent with 
the observed LRG power spectra, and also that the diagonal components of the covariance matrix from the mock 
catalogues give almost the same error as those of Eq. ((TE)) in the range of 0.02/iMpc-^ < ki< 0.2/iMpc"^ [39,,^]. 
Figure |6] shows the two dimensional map of the correlation matrix, 

,,(fc^, fc^.) ^ (26) 
y/ Cu{ki , ki)Cu{kj , kj ) 

for ^ = and 2 from 1000 mock catalogues. The binning of the covariance matrix is Afc = 0.01 ft,Mpc~^. One can see 
from Fig. IHlthat the off diagonal part is suppressed. 

The normalisation of the cosmological perturbations should be determined by the cosmic microwave background 
anisotropies, depending on the parameters n an kc of the f{R) model. However, the background expansion of the viable 
f{R) model is almost the same as that of the ACDM model, and the evolution of the matter density perturbations is 
only altered at late time, if compared with the ACDM model. This alteration will raise an additional integrated Sachs 
Wolfe effect on the CMB anisotropies due to the modified evolution of the matter density perturbation at late time. 
We neglect this effect on the normalisation of the perturbation, for simplicity. Then, we simply fixed the amplitude 
of the primordial cosmological perturbation by ag in the limit of large k^ i.e., the as of the ACDM model. 

Figure [S] shows that the shorter Compton wavelength model with a^ ~ 350km/s gives the best fit to the data. 
Figure [7] shows the contours of Ax^ on the kc — n plane. Here is computed with Eq. (p4)) with (|25l) by fitting the 
bias parameters and a^. The panels (a), (b), and (c) fix the normalisation of the perturbation to be as = 0.8, 0.82 
and 0.78, in the limit of large kc, respectively. The contour levels of Ax^ — 2.3 (dotted curve) and 6.2 (solid curve), 
correspond to la and 2a confidence, respectively. In figure [7] we used the Peacock and Dodds's formula (thin curve) 
and the Smith formula [82[ (thick curve), respectively, though the two curves almost overlap. The panels (a), (b), 
and (c) adopt the bias model of case 1. The panel (d) is the same as (a), but adopted the bias model of case 2. The 
left lower region in each panel is excluded. 

Figure [8] is the same as Fig. [71 but adopted the covariance matrix from the mock catalogues for the chi squared. 
The constraint of Fig. [5] is weaker compared with that of Fig. [7| Especially, the constraint for the model with larger 
n becomes weaker. However, Fig. [8] indicates that the long Compton wavelength case of the f{R) model with the 
smaller value of n is excluded. 

Thus far, we have used the redshift-space power spectrum ([T8l) . In order to check the reliability of our result, we 
next consider the other possible model for the redshift-space power spectrum, 

P^{k,fi) = {b\k)Pss{k) + 2fb{k)Ps0{k)fi^ + fPee{k)fi'') e-(-^''=^^^)', (27) 

where Pss{,k) is the nonlinear matter power spectrum, Peeik) is the power spectrum of the velocity divergence, and 
Pseik) is the cross power SMctrum of matter and the velocity divergence. This model is obtained from the model 
proposed by Scoccimarro [8q| and assumes a linear bias relation. Very recently, Jenning et al. proposed a fitting 
formula for the redshift-space power spectrum of the form (j27p . assuming b{k) = 1. The fitting formula relates the 
nonlinear matter power spectrum Pssik) to Pse{k) and Pgg{k). By using the N-body simulations it was demonstrated 
that the fitting formula is accurate to better than 10% for the ACDM model and quintessence dark energy models for 
k ^ 0.2/iMpc~^. Although the accuracy of the fitting formula for the /(i?) model has not been explicitly demonstrated, 
we assume its validity, and use it in the following Ax^ calculations. 

Figure [9] shows the contours of Ax^ on the kc — n plane, the same as Fig. [Tj but with the covariance matrix from 
the mock catalogues and the redshift space power spectrum (|27p . In the original formula, a^ is obtained from Pgg{k), 
however, we assumed a^ to be a fitting parameter, as is done in Fig. [5] This figure shows that the constraint becomes 
weaker when compared to the previous model (1181) . The models with large value of n are not constrained. However, 
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the long Compton wavelength case of the f{R) model with the smaller value of n is excluded. This new model predicts 
that Pse{k) is smaller than Pss{k) for values of fc O.lMpc"^, which reduces the quadrupole power spectrum and 
thus weakens the constraint. 

Let us compare our result with the other constraints on the f{R) model. Refs. [g^, [Ml have investigated the 
constraints on the f{R) model for the case n — 1/2. In Ref. [GOj], the constraint from the CMB anisotropics through 
the integrated Sachs Wolfe effect is investigated. However, the constraint is weak. Only the horizon-scale Compton 
wavelength model is excluded. In Ref. ^61'], the constraint from the cluster number count is investigated. Though it 
is restricted to the case n = 1/2, they obtained \fm\ ^ 10"'*, where fao is the value of fn at the present epoch. In 
the case n = 1/2, |/i^o| is related to fcc by 

/c, -0.04 -— hMpc-\ (28) 
\\JR0\J 

Ref. (HH reports that kc ^ 0.04ft,Mpc^^ is excluded. The constraint is similar to our result, when the redshift-space 
power spectrum ([T8|) is used (See Figure[8]). When arguably more accurate model p7| is used, the constraint becomes 
slightly weaker than that of ([TSl) (See Figure [9]). 



V. FUTURE PROSPECT OF MEASURING COMPTON SCALE 



In this section, we estimate future prospects of constraining the Compton scale with the use of the Fisher matrix 
technique, which is frequently used for estimating minimal attainable constraint on model parameters. We focus on 
the error of the Compton wavenumber kc- We adopt the Fisher matrix of the form (e.g., {87]). 



'"4^'A_ i-i de^ de^ (Pgai(fc,M) + i/")2' ^' 



where 9^ denotes a model parameter, V is a survey volume, n is a mean number density of galaxies. 

In the Fisher matrix analysis, for simplicity, we consider the 6 parameters kc, n, CTv, bo, hi and a, adopting the 
bias model of case 1. The panel (a) of Fig. 1101 shows the Icr error Afcc, in determining the Compton wavenumber 
kc as a function of the target value of fee, assuming a redshift survey like the SUMIRE (SUbaru Measurement of 
Imaging and REdshift of the universe) [9], which assumes the survey parameters like those of the WFMOS survey 
[ssj . the range of the redshift 0.9 < z < 1.6, the survey area 2000 square degrees, and the mean number density 
n = 4 X 10~'*(/i~iMpc)~^. Here we adopted the target values CTv — 400km/s, 6o = 2.5, hi — 0.5, a = 0.5, and n — 1/2, 
1, 2, and 4, from the bottom to the top, respectively. The other parameters are fixed fip = 0.28, h — 0.7, — 0.96, 
and the normalisation so as to be cts = 0.8 in the limit of the ACDM model. We obtained A/cc by marginalizing the 
Fisher matrix over the 5 parameters n, CTv, h^, hi and a. The panel (b) of Fig. 1101 shows the relative error Akc/kc. 

In the Fisher matrix we used the power spectrum in the range of wavenumbers k < 0.3/iMpc~^. This immediately 
implies that the redshift survey cannot be very sensitive to the models with the short Compton wavelength, as seen 
from figure [TOl The error becomes very large for kc ^ 0.1/iMpc^^, but it will be possible to obtain a useful constraint 
on the Compton scale, in principle, for models with kc ^ O.l/iMpc^^. However, the constraint becomes weak for the 
case of large n. 

The panel (a) assumes the power spectrum analysis without dividing the full galaxy sample, which spans the redshift 
range 0.9 < z < 1.6, into redshift subsamples. The panel (c) assumes the case when the galaxy sample is divided into 
the three subsample in redshift bins and that the power spectra are obtained from each subsample. In this case, the 
parameters a^, h^, hi and a should be fitted in each redshift bin, and the total number of parameters in the Fisher 
matrix analysis is 14. The panel (d) is the relative error, corresponding to (c). The cosmological parameters are 
the same as those of (a). The possible advantage of this method is that the additional information of the redshift 
evolution might improve the constraint. One can see that the constraint is improved in comparison with the panel 
(a) or (b). The degree of the improvement is small for n — 1/2, but is not negligible for the case n = 4. This is 
understood because the redshift evolution of the Compton scale is faster for larger n. 



VI. SUMMARY AND CONCLUSIONS 



In this paper, we determined a cosmological constraint on the viable f{R) model based on the redshift-space 
distortion by measuring the monopole and quadrupole spectra of the SDSS LRG sample of DR7. The monopole and 
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FIG. 10: (a) la error Afcc as a function of the target value of fee. The result is based on the Fisher matrix analysis with the 6 
parameters, fee, n, (Jv, and bo, foi and a for the bias model 1, and marginalized over the 5 parameters other than kc- The target 
parameters are 6o = 2.5, b\ — 0.5, a = 1/2, and n is chosen n = 1/2, 1, 2, 4 from the bottom to the top, respectively. The 
other parameters are fixed Qq — 0.28, h — 0.7, Us = 0.96, and the normalisation erg — 0.8 in the limit of the ACDM model. 
Eq. (|18p with the Peacock and Dodds nonlinear fitting formula is adopted, (b) the relative error Afcc/fcc corresponding to (a), 
(c) and (d) are the same as (a) and (b), respectively, but assumed the analysis where the full sample is divided into 3 redshift 
bins. 

the quadrupole spectra are used to fit the bias parameters and to constrain the growth factor and the growth rate of 
the density perturbations, which depend on the Compton scale of the f{R) model. 

Our results show that short Compton wavelength model fits the data better, while the long Compton wavelength 
model is excluded, though the constraint depends on the evolution parameter n. For the case n = 1/2, our constraint 
is similar to that from the cluster number counts reported in [61[. When we adopt more accurate model for the 
redshift-space power spectrum [s^l, the constraint becomes slightly weaker. However, the long Compton wavelength 
case of the f{R) model with the smaller value of n is excluded. Our results exemplify that the redshift-space distortion 
is quite useful in testing gravity theory. We also demonstrated that a future redshift survey like the WFMOS/SUMIRE 
is potentially useful in obtaining a constraint on the Compton wavelength scale. 

We acknowledge that the widely used theoretical model of the anisotropic power spectrum adopted in the present 
paper might need careful improvements. We adopted the Peacock and Dodds formula and the Smith formula for the 
nonlinear modelling of the mass power spectrum. Our results do not significantly depend on the choice. However, there 
might be a need to adopt a more sophisticated formula for the precise nonlinear modelling within the framework of 
the modified gravity, as demonstrated by Koyama, Taruya, Hiramatsu [11] . The treatment of the Finger of God effect 
in our paper was simple, which assumed the exponential distribution function for the pairwise velocity and introduced 
one free parameter - the pairwise velocity dispersion. In reality it might not be an adequate model to describe the 
nonlinear region of the redshift-space power spectrum [85|. We checked the reliability of our results by adopting the 
other possible model proposed in Ref.[84|, extensively applying the fitting formula to the f{R) model, whose accuracy 
in this case, however, has not been demonstrated. We found that there is a non-negligible effect on the constraint on 
the f{R) model. Therefore, a more precise modelling of the redshift-space power spectrum should arguably be needed 
in the future. Concerning the modelling of the clustering bias, we adopted a simple scale-dependent bias. Here too 
there is potentially a lot of room for improvement. These issues are out of scope for the present paper, but need to 
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be elaborated for a precise test of gravity with the future redshift surveys. 
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